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Abstract 

We consider the neutrino oscillations within the minimal supersymmetric standard 
model with R-parity violation (^f p MSSM). The Super-Kamiokande atmospheric neu¬ 
trino data are used to set limits on the bilinear R-parity violating terms of the ^? p MSSM. 
These very stringent limits are out of reach of the other experiments at present and in 
the near future. 
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The neutrino masses and mixing constitute one of the most pressing problem of particle 
physics. At present, there is a convincing experimental evidence for a non-trivial structure 
of the three-generation neutrino mass matrix. It derives from the solar [l| and atmospheric 
|| neutrino deficit as well as, likely, from the LSND experiment ||. All the three types 
of observations can be explained by the (anti)neutrino oscillations. An experimental break¬ 
through in exploration of the atmospheric neutrino anomaly has been recently made by 
the Super-Kamiokande collaboration [Q measured the zenith angle dependence of the atmo¬ 
spheric neutrino flux. This high statistics experiment confirmed that the most likely solution 
of the atmospheric neutrino anomaly implies neutrino oscillation preferably in z/ p —* u T chan¬ 
nel ||, |5| or, maybe, partially in active-sterile channel —► u s ||. This result may have far 
reaching consequences for the new physics phenomenology j7]. 

From the theoretical side the neutrino mass problem has been attracting many efforts 
during all the history of modern particle physics. It stimulated construction of many inter¬ 
esting approaches which have their meaning far beyond the study of the neutrino proper¬ 
ties. The presently popular supersymmetric (SUSY) models with R-parity non-conservation 
(Ifl p SUSY) may provide us with a new insight into the neutrino mass problem. In this 
framework neutrinos can acquire the tree-level supersymmetric masses via the mixing with 
the gauginos and higgsinos at the weak-scale || |DJ which is possible due to the presence 
of the bilinear R-parity violating operators in the superpotential and soft SUSY breaking 
sector. Remarkable, that this mechanism does not involve the physics at the large energy 
scales Mi nt , ~ 0( 10 12 GeV), in contrast to the see-saw mechanism, but relates the neutrino 
properties to the weak-scale physics accessible for the experimental studies. 

In the present paper we consider an impact of the new atmospheric neutrino data on 
the bilinear R-parity violation. We use the constraints on neutrino oscillation parameters 
recently obtained in Ref. J5|] from the combined fit to the Super-Kamiokande atmospheric 

On this basis we derive new stringent 


and to the CHOOZ reactor neutrino data 10 


constraints on the R-parity violating parameters of all the three generations. These con¬ 
straints may have important implications for the lepton number/flavor violating processes 
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such as the neutrinoless double beta (Ov/3/3) decay, muon-electron (positron) conversion, etc. 
The djt p SUSY contribution to these processes can be expressed in terms of the quantities 
observed in the neutrino oscillation experiments. We have found that the Super-Kamiokande 
casts significantly more stringent limits on the 1st generation bilinear R-parity violation than 
the present Ozz/3/3-decay experiments. Previously, before the Super-Kamiokande data have 
been published, Oz//?/?-decay experiments provided the best limits of this type 
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Let us consider the neutrino oscillations in the minimal supersymmetric standard model 
with R-parity violation (^? p MSSM). Recall that the R-parity is a multiplicative Z 2 symmetry 


3B+L+2S 


where S, B and L are the spin, the baryon and the lepton 


defined as R p = (—1) 
quantum numbers. 

The R p violation is introduced into the theory through the superpotential 

WV = hjkULjEi + KjkUQjDl + ^L S H 2 + \" jk U?D c j D e k , (1) 

and soft SUSY breaking scalar potential terms 

VS* = A ijtUijEi + K, t UQM + A " iik vtb c M + &1,h 2 + p' Lj H\ + h.c. ( 2 ) 


Here L , Q stand for lepton and quark doublet left-handed superfields while E c , U c , D c for 
lepton and up, down quark singlet superfields; H\ and Ho are the Higgs doublet superfields 
with a weak hypercharge Y = — 1 , + 1 , respectively. Fields in Eq. (||D are the scalar 
components of these superfields. Summation over the generations is implied. To prevent 
fast proton decay it is a common practice to set A" = A" = 0. 

The impact of the R-parity violation on the low energy phenomenology is twofold. First, 
it leads to the lepton number/flavor violating interactions directly from the trilinear terms in 
Wfi p . Second, bilinear terms in W^ p and in V^° ft generate the non-zero vacuum expectation 
value for the sneutrino fields (zJj) 7 ^ 0 and cause neutrino-neutralino as well as electron- 
chargino mixing. The mixing gives rise to the non-trivial neutrino mass matrix and brings 
in the new lepton number/flavor violating interactions in the physical mass eigenstate basis. 
These interactions contribute to some exotic processes such as Ov/3/3-decay, — e ± con¬ 


version, etc. JTI| and, therefore, can be constrained experimentally. In the present paper 
we analyze constraints on the neutrino mass matrix of the dji p MSSM from the neutrino 
oscillation experiments. 

Let us recall the mechanism of neutrino mass generation in the Jfip MSSM. The bilinear 
terms in the Eqs. (|l|) and (|2j) lead to the terms in the scalar potential linear in the sneutrino 
fields As a result, at the minimum of the potential these fields acquire non-zero vacuum 
expectation values ( i>i) 7 ^ 0. Thus, the MSSM vertices Zvv and Wev create the gaugino- 
lepton mixing mass terms Zv(v), We(v) (with IU, Z being wino and zino fields). Mixing of 
the leptons with the Higgsinos comes from the /ijLjLU term in the superpotential ([[]). 

In the two component Weyl basis 


vU T 


(0) V'e) " \i 

the mass term of the neutral fermions is 

1 


= ( 1 //, v' T , -iA, -iA 3 , H\\ Hi), 


= -f^o^o) - H.c., 


( 3 ) 


( 4 ) 
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Here v' e ,i/ ,j/ t are the neutrino weak eigenstate fields, A and A 3 are the Lfiy and SU 2 l 
gauginos, respectively, while the higgsinos are denoted as H® 2 . 

The 7x7 mass matrix has the see-saw structure 


M 0 = 



m 

M x 


with 3x4 matrix 


/ 

m = 

V 


—M z s w cpv e 
-M z s w cpv M 
—M z s w cpv r 


M z c w cpv e 

M Z C w CgV fl 

M z C w CgV r 


0 -/m e \ 

0 -f.IU M , 

0 — IXU T ) 


(5) 


( 6 ) 


originating from the ^bilinear terms in the superpotential (jTj) and the soft SUSY breaking 
sector (^). Here tan/? = (H®) / (H °), sw = sin0 w , cw = cos#^, sp = sin/?, cp = cos/?. 
The supersymmetric Higgs mass /x occurs in ^lHiH 2 term of the MSSM superpotential [[12] . 
We introduced the quantities v a = (£> a )/(i7°), u a = /i a /n which are assumed to be small, 
v a ,u a « 1, and used as expansion parameters. For the details of this approximate approach 
we refer to Ref. [[□[]. Hereafter a = e, /i, r is a flavor index. In Eq. 

MSSM neutralino mass matrix in the basis {— iX, — xA 3 , H®, H®} 


M x is the usual 4x4 


ilW = 


^ Ml 0 —M z S\yCp M z S\ySp ^ 

0 M 2 M z C\yCp —M Z C W Sp 

—M z swCp M z cwCp 0 —/x 

V M Z S\ySp —M Z CwSp —yU 0 ) 


(7) 


In the mass eigenstate basis defined as 


*(0)i = 3y*'( 0)i , 

the 7x7 neutral fermion mass matrix Ado in Eq. (^|) becomes diagonal 

5*7Vd 0 S t = Diagjm^, }, 


( 8 ) 

(9) 


where m Ui and m Xi are the physical neutrino and neutralino masses. For the considered case 
of the tree-level mass matrix only one neutrino has a non-zero mass m Vl = m„ 2 = 0, m U3 ^ 0. 
This is a result of the minimal field content and the gauge invariance. The mass matrix A4 0 
in Eq. (|5]) has such a texture that its first three rows and the last one are linearly dependent 
and, therefore, two neutrino mass eigenstates are degenerate massless states. The v\ — v 2 
mass degeneracy is lifted by the 1-loop corrections as well as by the non-renormalizable terms 
in the superpotential. As a result z/ 1)2 acquire small non-equal masses m Ul ^ m V2 ^ 0 |J. 
This effect is not important for our analysis since we confine ourselves to the atmospheric 
neutrino anomaly explained as —> v T oscillations. In this case the 1-loop corrections and 
non-renormalizable terms can be neglected as they give only small corrections to the tree 
level value of the corresponding oscillation parameter Srri 2 tm = m 2 ui 2 — m 2 3 = m 2 3 ^ 0. This 
is not the case for the treatment of the solar neutrino data, interpreted in terms of v e —> z/ /t 
oscillations with Sm 2 ol = m 2 — rri 2 2 = 0 at tree level in the Ifl p MSSM. 
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The mixing within the neutrino sector to leading order in the small expansion parameters 
v a ,u a is described by 

= V£ r i/ a ( 10 ) 

with v' a = (i/, 1 / , v' T ) being the weak basis. 

The 3x3 matrix diagonalizes the ft p -induced effective neutrino mass matrix: 

V^ T m eS yM = Diag{0, 0, m U3 }. (11) 

The tree-level expression for this mass matrix can be found in Ref. 0 - The only non-zero 
eigenvalue which we identify with the tau neutrino mass is given by 

~ ^ IA 12 


m. 


' V3 -‘mssm 


■EIA. 

a 


where 


The function F MSSM depends only on t. 
rninant of the MSSM neutralino mass 


4 Det M x 

MSSM parameters. 


a _ \ , 1 j\.. v _ „ 2 M ! + tan 2 9 w M 2 

A-a 9\^cx) \Hi)fl a ^ -^MSSM 92 ^ ■ - - 

depends only on the R-parity conserving 
vf neutralino mass matrix (|7|) is 

Det M x = sin 2 /3 //(Mi + tan 2 9 W M 2 ) — MiM 2 /i 2 . 

The ft p MSSM neutrino mixing matrix can be written in the following form [] 

0 


( 12 ) 

(13) 
. The deter- 

(14) 


yM 


(15) 


/ cos 0 i 3 0 — sin 0 13 

sin 023 sin 0 13 cos 0 2 3 sin 0 23 cos 0 i3 
\ cos 0 23 sin 013 — sin 0 23 cos0i 3 cos0 23 

where the mixing angles are defined as follows: 

tan 0i 3 =- = , tan0 23 = — J - . 

yyTAi' a t 

We used a freedom in definition of the mass eigenstates v \, with m ui 
V$ = 0 in Eq. ®D. 

Now, let us turn to neutrino oscillations in vacuum. The survival probability for a given 
neutrino flavor v a is described by the formula [|b| 

„ „ 5m% ■ L 

= 1 - 4 Z PtxjPcA sm , 

i<j 

r(y) 
v aj 


(16) 

= m U2 = 0 to set 


P(i 


v a -> v a 


( 17 ) 


where Sm 2 j = ml. — ml. and P a j = V aj 

lr The ^pMSSM neutrino mixing matrix given in Ref. 113 contains misprints. 
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From the combined fit to the Super-Kamiokande jj and the CHOOZ |10] data the fol¬ 
lowing constraints at 95% C.L. were recently found in Ref. ||: 


0.5 

< 

<5™atm/( 10 ~ 3 eV 2 ) < 10.0, 

(18) 

0.00 

< 

Pe3< 0.08, 

(19) 

0.25 

< 

P^3 <0.75, 

(20) 

0.25 

< 

Pt3 < 0.75. 

(21) 


The analysis of Ref. |J was made under the assumption |<5m 2 tm | » |<5m 2 ol |, with <5m 2 tm = 
5m\ 2 « and Sml ol = This covers a particular case of the ]Jl p MSSM neutrino 

mass spectrum Eqs. (|TT|) and (|l2|) which leads to the relations 


( 22 ) 


= 5m 32 = = ™v 3 = F MSS m 



and 5rnl tm » Sm 2 sol = m 2 2 


m 2 since m Vl 2 are zero at the tree level. 


The mixing parameters in Eqs. 
)|) as 


are determined by the tree level mixing matrix 


Pa3 ~ 


A 2 


Ea|A Q | 2 ' 

It is useful to introduce the effective SUSY parameters 

A — A 2 . tp — p 

'/'■Cv x 1 ^v/ - 1 - MQCM 1 ( 


«3 


dm 2 


atm • 


(23) 


(24) 


An advantage of these parameters for our analysis follows from the fact that they depend 
on the ]Jl p parameters (z/ a ), /i Q from only one generation a unlike parameters P a3 , 5m 2 atm , 
depending on the ]fi p parameters from different generations. 

The Eqs. (P~8|)—(|2T|) cast the constraints on the effective parameters A a : 


o < Ae < 8 • 10" 3 eV, 


5.6 • 10“ 3 eV < A„, T < 7.5 


KT 2 eV. 


(25) 


This constraint represents a complex exclusion condition imposed by the neutrino oscillation 
data on the lji p MSSM parameter space. 

It is instructive to obtain typical individual constraints on the superpotential parameters 
fj, a and on the sneutrino vacuum expectation value (P a ). These constraints can be derived 
from Eqs. (P5| ) at typical weak scale values of the MSSM parameters M li2 ~ /i ~ M SVSY and 
tan/? = (H 2 ) / (Hi) = 1. Here M SVSY is a characteristic SUSY breaking mass scale varying 
most likely in the interval 100 GeV $ M SVSY < 1 TeV motivated by non-observation of the 
SUSY particles and by the ’’naturalness” arguments. Following common practice, we also 
assume absence of a significant cancellation between the two terms in definition of A a given 
in Eq. (|I3]). Thus, we come up with the following constraints 


I he | > |(^ e )| 


(A4usy/ 100 GeV) 


< 90 KeV, 76 KeV < 


I AV,t| ; | (G^t) I 

(A4usy/ 100 GeV) 


< 276 KeV. 


(26) 
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These constraints apply for all the quantities /r a , (n Q ) separated by commas in the above 
formulas. To our knowledge these stringent limits have not been previously derived in the 
literature. One can find in the published papers only those constraints which involve the 
combinations of the 1st and the 2nd generation bilinear Jfl v parameters |T3| or contain only 
the 3rd generation ones [HJ. The constraints from Refs. |H|, JT5[ are significantly weaker 


(by several orders) than those in Eq. (^6|). Recent study in Ref. [16] concentrate mostly on 


a possibility of simultaneous solution of all the three known neutrino anomalies in the GUT 
constrained lfi p MSSM scenarios rather than on extraction of constraints on the low energy 
parameters like in Eqs. ( |25|) and (|25|). 

It is interesting to note that in the lfi p MSSM the generation average Majorana neutrino 
mass is given at tree level by the formula: 


(™>v) = J2 m vi ( V e 


ei 


= P p *\ 5ml 


6 atm 


(27) 


Thus, it is expressed in terms of the neutrino oscillation parameters. 

As is known the parameter {m v ) plays a clue role in the phenomenology of the Ozyd/3-decay. 
It determines the classical Majorana neutrino exchange contribution to the Ozyd/3-decay half- 
life as 




0+r 1 = G m \M. 


to 

K) 


m e 


(28) 


Here Goi and M. v are the phase space factor and the nuclear matrix element. For O^/3/3-decay 
of 73 ' 


j Ge we have G 01 = 7.93 ■ 10" 15 years" 1 


and M. v = —3.4 11 


The nuclear matrix ele¬ 
ment is calculated in the pn-RQRPA approach [|T8|. The contribution ([28]) represents a dom¬ 
inant effect of the bilinear R-parity violation in 0z//l/3-decay |TT 
limit on Ozyd/3-decay half-life of 76 Ge 


g T 1 ° / f /3 (o+ o + ) > 


1/2 


From the best experimental 
1.1 x 10 25 years (90% C.L.) 


one obtains {m v ) < 0.50 eV. From the Super-Kamiokande and CHOOZ constraints in Eq. 
( :m-m we derive significantly more stringent limit 


< 0.80 


10 -2 eV 


(29) 


Thus we conclude that at present the neutrino oscillation experiments are more sensitive to 
the bilinear R-parity violation than the 0uf3 /9-decay experiments. This conclusion holds for 
the whole parameter space of the MSSM since both type of the experiments place constraints 
on the same effective parameter A e which absorbs the dependence on the MSSM parameters. 
To be comparable in sensitivity to this parameter with the present neutrino oscillation exper¬ 
iments the O/yd/i-experinients have to reach the half-life lower limit ~ 10 29 years. This 

is unrealistic for the running experiments and might only be possible for recently proposed 
large scale germanium detector 0. 
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